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Low-Energy Supersymmetry and the Tevatron Bottom-Quark Cross Section 
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A long-standing discrepancy between the bottom-quark production cross section and predictions of 
perturbative quantum chromodynamics is addressed. We show that pair production of light gluinos, 
of mass 12 to 16 GeV, with two-body decays into bottom quarks and light bottom squarks, yields a 
bottom-quark production rate in agreement with hadron collider data. We examine constraints on 
this scenario from low-energy data and make predictions that may be tested at the next run of the 
Fermilab Tevatron Collider. 
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The measured cross section for bottom-quark produc- 
tion at hadron collider energies exceeds the expecta- 
tions of next-to-leading order calculations in perturba- 
tive quantum chromodynamics (QCD) by about a factor 
of 2 The next-to- leading order (NLO) corrections 
are large, and it is not excluded that higher order effects 
in production or fragmentation may resolve the discrep- 
ancy. However, this long-standing discrepancy has so far 
resisted satisfactory resolution within the standard model 
(SM) of particle physics g]. The disagreement is surpris- 
ing because the relatively large mass of the bottom quark 
sets a hard scattering scale at which perturbative QCD 
computations of other processes are generally successful. 
In this Letter we explore an explanation of the discrep- 
ancy within the context of the minimal supersymmetric 
standard model (MSSM) Q . We postulate the existence 
of a relatively light gluino g (mass ~ 12-16 GeV) that 
decays into a bottom quark and a light bottom squark b 
(mass ~ 2-5.5 GeV). The g and the b are the spin-1/2 
and spin-0 supersymmetric partners of the gluon (g) and 
bottom quark (6). In our scenario the b is either long- 
lived or decays hadronically. We obtain good agreement 
with hadron collider rates of bottom-quark production. 
Several predictions are made that can be tested readily 
with forthcoming data from run II of the Fermilab Teva- 
tron collider. 

Our assumptions are consistent with all experimental 
constraints on the masses and couplings of supersymmet- 
ric particles PHl^]- An analysis of four-jet data by the 
ALEPH Collaboration disfavors a g with mass rrig < 6.3 
GeV [(§ but does not cover the mass range considered in 
this Letter. An analysis by the UA1 Collaboration || 
of the mass range 4 < nig < 53 GeV does not apply to 
our scenario because they assume that the gluino decays 
into two quarks plus large missing energy. However, a 
new comparison to the UA1 data [Q does see an excess 
in the 6-quark cross section, as would be expected in our 
model. 



If the light b is an appropriate admixture of left-handed 
and right-handed bottom squarks, its tree-level coupling 
to the neutral gauge boson Z can be small, leading to 
good agreement with the Z-peak observables ||. Bot- 
tom squarks make a tiny contribution to the inclusive 
cross section for e + e~ — * hadrons, in comparison to the 
contributions from quark production, and bb resonances 
are likely to be impossible to extract from backgrounds 
. One can study the angular distribution of hadronic 
jets produced in e + e~ annihilation in order to bound 
the contribution of scalar-quark production. Spin-1/2 
quarks and spin-0 squarks emerge with different distri- 
butions, (1 ± cos 2 #), respectively. We refit the angular 
distribution measured by the CELLO Collaboration || 
and find it is consistent with the production of a sin- 
gle pair of charge-1/3 squarks along with five flavors of 
quark- antiquark pairs. The exclusion by the CLEO Col- 
laboration jl0| of a b with mass 3.5 to 4.5 GeV does not 
apply since that analysis focuses on the decays 6 — * civ 
and b — > cl [ pd| . Thus, measurements at e + e~ colliders 
do not significantly constrain b masses in the region of 
interest. 

A long-lived b is not excluded by conventional searches 
at hadron and lepton colliders, but an analysis similar to 
that of Ref . Jl2j should be done to verify that there are no 
additional constraints on the allowed range of b masses 
and lifetimes. Alternately, the b could decay hadronically 
via a baryon-number- and i?-parity-violating interaction 
into soft light hadrons which will fall within the cone of 
the b jet. The main constraint is that the b does not 
decay into hard leptons or leave a large missing energy. 
The b and g masses we consider are also compatible with 
theoretical constraints which require the absence of color 
and charge breaking minima in the scalar potential ]13|| , 

Because the excess production rate is observed in all 
bottom-quark decay channels, an explanation in terms of 
"new physics" is guided towards hypothesized new par- 
ticles that decay either like bottom quarks or directly 
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to bottom quarks. The former is difficult to implement 
successfully in the MSSM, as explained in the discus- 
sion of alternative scenarios at the end of the Letter. 
We adopt the latter. In our scenario, light gluinos are 
produced in pairs via standard QCD subprocesses, dom- 
inantly g + g — > g + g at Tevatron energies. The g has 
a strong color coupling to 6's and 6's and, as long as its 
mass satisfies nig > nib + mg, the g decays promptly to 
6 + 6. The magnitude of the 6 cross section, the shape 
of the 6's transverse momentum pxb distribution, and 
the CDF measurement Oj of B° — B a mixing are three 
features of the data that help to establish the preferred 
masses of the g and 6. 

Including contributions from both q + q — > g + g and 
g + g — > g + g |l5|, we show in Fig. 1 the integrated pxb 
distribution of the 6 quarks that results from g — > 6 + 6, 
for m g =14 GeV and mg = 3.5 GeV. The results are 
compared with the cross section obtained from NLO per- 
turbative QCD Q a nd CTEQ4M parton distribution 
functions (PDF's) |7|, with m b = 4.75 GeV and a re- 
normalization and factorization scale /i = \/m\ + p\ h - 
SUSY-QCD corrections to 66 production are not included 
as they are not available and are generally expected to 
be somewhat smaller than the standard QCD corrections. 
A fully differential NLO calculation of g-pair production 
and decay does not exist either. Therefore, we compute 
the g-pair cross section from the leading order (LO) ma- 
trix element with NLO PDF's 0, ^ = «/m| + a 

two- loop a g , and multiply by 1.9, the ratio of inclusive 
NLO to LO cross sections Q. 

A relatively light gluino is necessary in order to obtain 
a bottom-quark cross section comparable in magnitude to 
the pure QCD component. Values of m g ~ 12-16 GeV 
are chosen because the resulting g decays produce prb 
spectra that are enhanced primarily in the neighborhood 
of Px™ — m g where the data show the most prominent 
enhancement above the QCD expectation. Larger values 
of frig yield too little cross section to be of interest, and 
smaller values produce more cross section than seems tol- 
erated by the ratio of like-sign to opposite-sign leptons 
from 6 decay, as discussed below. The choice of mg has 
an impact on the kinematics of the 6. After selections 
on Py™, large values of mg reduce the cross section and, 
in addition, lead to shapes of the pTb distribution that 
agree less well with the data. The values of mg and nig 
are correlated; similar results to those shown in Fig. 1 
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can be obtained with nig 



12 GeV, but nil — m t>- 



After the contributions of the NLO QCD and SUSY 
components are added (solid curve in Fig. 1), the mag- 
nitude of the bottom-quark cross section and the shape 
of the integrated p™ l b n distribution are described well. A 
theoretical uncertainty of roughly ±30% may be assigned 
to the final solid curve, associated with variation of the 
6 mass, the scale, and the parton distributions. 
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FIG. 1. Bottom-quark cross section in pp collisions at 
■\/S — 1.8 TeV for prb > Prt" with a gluino of mass m,g — 14 
GeV and a bottom squark of mass = 3.5 GeV. The dashed 
curve is the central value of the NLO QCD prediction. The 
dotted curve shows the pr spectrum of the b from the super- 
symmetry (SUSY) processes. The solid curve is the sum of 
the QCD and SUSY components. Data are from Ref. 

The SUSY process produces bottom quarks in a four- 
body final state and thus their momentum correlations 
are different from those of QCD. Angular correlations 
between muons that arise from decays of 6's have been 
measured Jl4| , ^9| . Examining the angular correlations 
between 6's in the SUSY case we find they are nearly 
indistinguishable from those of QCD once experimental 
cuts are applied. 

Since the g is a Majorana particle, its decay yields 
both quarks and antiquarks. Gluino pair production and 
subsequent decay to 6's will generate 66 and 66 pairs, as 
well as the 66 final states that appear in QCD production. 

We perform an exact matrix-element calculation of 
the four-body cross section for like-sign and opposite- 
sign bottom quarks from g pair-production and decay. 
When a gluino is highly relativistic, its helicity is nearly 
the same as its chirality. Therefore, selection of g's 
whose transverse momentum is greater than their mass 
will reduce the number of like-sign 6's. In the limit 
of either massless g's or very high pt cuts, the like- 
sign to opposite-sign ratio reduces to y/(l — y) where 
V = \ sin 2 26>g and sin 9^ is the left-handed component of 
the lightest bottom-squark mass eigenstate. There is a 
strong suppression of like-sign 6's if the mixing angle is 
small. For the case under consideration, the mixing of the 
bottom squark is determined by the condition that the 6 
coupling to the Z boson is small Q , namely, sin 6>g — 0.38. 
In the intermediate pr region, however, the like-sign sup- 
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pression is reduced. The cuts chosen in hadron collider 
experiments for measurement of the like-sign to opposite- 
sign muon ratio result in primarily unpolarized g's, and, 
independent of the 6 mixing angle, an equal number of 
like-sign and opposite-sign 6's is expected at production. 

The SUSY mechanism leads to an increase of like-sign 
lcptons in the final state after semi-leptonic decays of 
the b and b quarks. This increase could be confused with 
an enhanced rate of B° — B° mixing. Time-integrated 
mixing analyses of lepton pairs observed at hadron col- 
liders determine the quantity x, expressed conventionally 
as X = fdXd + fsXs, where f d and f s are the fractions of 
B® and B® hadrons, respectively, in the sample of semi- 
leptonic B decays, and X/ is the time- integrated mixing 
probability for B®. The quantity 2x(l — X) is the fraction 
of 66 pairs that decay as like-sign b's. Our SUSY mech- 
anism can be incorporated by introducing Xcff such that 
2Xcff(l - Xcff) = [2x(l ~X) + G/2]/(l + G), where G is 
the ratio of SUSY and QCD bottom-quark cross sections 
after cuts. The effective mixing parameter constrains G: 



Xcff 



X 



Vl + G 



1 - 



1 



VT+G 



(1) 



To estimate X e ff> we assume that the world average 
value x — 0.118 ± 0.005 |g0[ represents the contribution 
from only the pure QCD component. We determine the 
ratio G in the region of phase space where the measure- 
ment is made |Q, with both final b's having pt of at 
least 6.5 GeV and rapidity \yb\ < 1. The ratio is com- 
puted with LO matrix elements, NLO PDF's, a s at two 
loops, and fi — m x , where x is b or g. The SUSY and 
QCD cross sections are multiplied by 1.9 and 2.3, re- 
spectively, to account for NLO effects llqJlql . For gluino 



masses of 



= 14 and 16 GcV, we obtain G = 0.37 



and 0.28, respectively, with = 3.5 GeV. We compute 
Xcff — 0.17 for rrig — 14 GeV, and XeB — 0-16 with rrig = 
16 GeV. 

To estimate the uncertainty on G, we vary the scale at 
which the cross sections are evaluated between fi = m x /2 
and ii = 2m x . Uncertainties of ±50% are obtained and 
lead to uncertainties in the determined values of 5xeS — 
±0.02. Additional uncertainties arise because there is no 
fully differential NLO calculation of gluino production 
and subsequent decay to b's. 

Our expectations may be compared with the CDF 
Collaboration's published value Xcff = 0.131 ± 0.02 ± 
0.016 @. Values of m g > 12 GeV lead to a calculated 
Xeff that is consistent with the measured value within 
experimental and theoretical uncertainties. 

The effects of light gluinos and bottom squarks on the 
strong gauge coupling a s are potentially very significant. 
In the SM, a global fit to all observables provides an indi- 
rect measurement of a s at the scale of the Z boson mass 
M z . The value a s (M z ) ~ 0.119 ± 0.006 describes most 
observables properly |2M. QCD induced processes such 



as jet cross sections at hadron colliders and the top-quark 
cross section are consistent with this value of a s (Mz). A 
light g with mass about 15 GeV and a light b modify 
a s (Mz), determined by extrapolation from experiments 
performed at energies lower than m g . The result is a 
shift of 0.007 in a s derived from these experiments, to 
a s (M z ) = 0.125. The value of a s (M z ) obtained from 
the Z hadronic width is modified by the presence of light 
6 s and light g's, but the effect is slight because these su- 
perpartners approximately decouple from the Z. In the 
presence of light g's, all values of a s (Mz) still fall within 
the range of experimental uncertainty, with a slight pref- 
erence for the upper edge of the range. We do not at- 
tempt a global fit to the data. 

Small values of rrig and of the masses of the third gener- 
ation squarks are helpful in supersymmetric solutions of 
the hierarchy problem of the SM. They reduce the fine- 
tuning needed, and avoid a color-breaking vacuum, in 
theories where SUSY is broken at the Planck scale Jl3|,^2| . 
Light gluinos are also helpful in improving gauge coupling 
unification by providing a light threshold in the evolution 
of a s [^3| . For gluinos of mass several hundred GeV, the 
strong coupling at Mz predicted from unification is some- 
what above 0.13. However, for gluinos of mass 15 GeV, 
this prediction becomes a s (Mz) ^ 0.127, in agreement 
with the measured value. The upper limit is saturated 
when the masses of the weak gauginos are of the order of 
the weak scale. 

Among the predictions of this SUSY scenario, the most 
clearcut is pair production of like-sign charged B mesons 
at the Tevatron collider. A very precise measurement 
of x m run II is obviously desirable. Since the fraction 
of 6's from gluinos changes with pTb, we also expect a 
change of x with the cut on pxb- Possible bound states 
of bottom squark pairs are not expected to be prominent 
in e + e~ annihilation but could be seen as mesonic 
resonances in 77 reactions or as narrow states in /i pair 
production in hadron collisions. The existence of light 6's 
means that they will be pair produced in partonic pro- 
cesses, leading to a slight increase in the hadronic dijet 
rate. Our approach increases the 6 production rate at 
DESY ep collider HERA and in photon-photon collisions 
at CERN Large Electron-Positron Collider (LEP) by a 
small amount, not enough perhaps if early experimen- 
tal indications in these cases are confirmed ^4j,^5|, but 
a full NLO study should be undertaken. Finally, to sat- 
isfy constraints from electroweak measurements, a light 
superpartner of the top quark, with mass about the top- 
quark mass should be present in the spectrum The 
lightest Higgs boson should be lighter than 123 GeV, and 
it should decay mainly into light 6's. A study of possible 
6 decay modes could be undertaken to see if our scenario 
can be made consistent with the Higgs-candidate events 
at LEP @. 

In our investigation, we consider and discard alterna- 
tive scenarios. If only the 6 is light, with 6 decay products 
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similar enough to those of the b quark, one can imagine 
that a light b might be sufficient for a good description 
of the Tevatron b cross section. The cross section for 
pair production of light &'s becomes comparable to the 
bottom-quark rate for m~ h ~ 3 GeV. A b with mass of 
about 3 GeV decaying through an i?-parity violating pro- 
cess into a r lepton and a charm quark could lead to an 
excess of the apparent b cross section. This scheme proves 
difficult to implement. Such a light b cannot decay to a 
J I ip and will therefore not reproduce the observed excess 
b cross section in this decay channel. Moreover, light fe's 
do not reproduce the observed pt dependence of the b 
cross section. 

If mj + mi > rrig, but m& < mg < m§, the gluino can 
decay into a b and a strange or down quark, followed by 
decay of the b into a b and a light neutralino x°. For 
this scenario, there must be a flavor violating coupling 
of g — b — s of about to suppress the branching ra- 
tio of the alternative g decay channel: g — ► gx°- The 
problem is that a light b proceeding from g decay, and 
decaying into a b and missing energy, would lead to a bb 
plus missing energy cross section of order of 100 pb, so 
large that it should have been observed in the data from 
run I of the Tevatron collider 27 1. 

In this Letter, we propose an interpretation of the 
excess bottom-quark production rate at the Tevatron 
that involves new physics and leads to several testable 
consequences. We show that pair production of low- 
mass gluinos that decay into bottom quarks and bot- 
tom squarks provides a good description of the b cross 
section and reproduces the measured ratio of like-sign 
to opposite-sign leptons. Our assumptions are consis- 
tent with all experimental constraints on SUSY particle 
masses. 
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